Dipole responses in Nd and Sm isotopes with shape transitions 
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Photoabsorption cross sections of Nd and Sm isotopes from spherical to deformed even nuclei are 
systematically investigated by means of the quasiparticle-random-phase approximation based on the 
Hartree-Fock-Bogoliubov ground states (HFB+QRPA) using the Skyrme energy density functional. 
The gradual onset of deformation in the ground states as increasing the neutron number leads to 
characteristic features of the shape phase transition. The calculation well reproduce the isotopic 
dependence of broadening and emergence of a double-peak structure in the cross sections without 
any adjustable parameter. We also find that the deformation plays a significant role for low-energy 
dipole strengths. The El strengths are fragmented and considerably lowered in energy. The summed 
El strength up to 10 MeV is enhanced by a factor of five or more. 
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Density functional theory has been widely used to de- 
scribe a variety of quantum many-body systems [ij in- 
cluding nuclear many-body systems Q • Recent advances 
in the computing capability together with the highly- 
developed techniques in the nuclear energy-density- 
functional (EDF) method allow ns to calculate the 
ground-state properties of nuclei in the entire mass re- 
gion Q. The nuclear ground-state deformation is one of 
them, which is an example of the spontaneous break- 
ing of rotational symmetry. Experimental evidences 
of the nuclear shape changes are related to low-lying 
quadrupole collectivity, such as the ratio of the excitation 
energy of 2+ and 4+ states i?4+ / E2+ , the reduced tran- 
sition probability B{E2\2^ O'^), etc. However, it is 
known that the nuclear deformation also affects the high- 
frequency collective modes of excitation, giant resonances 
(GRs) [ii. For instance, the peak splitting of the giant 
dipole resonance (GDR), which is caused by the differ- 
ent frequencies of oscillation along the long and short 
axes, has been observed in experiments Q . A typical ex- 
ample of the shape phase transition form spherical to de- 
formed ground states and an evolution of the deformation 
splitting in GDRs have been observed in Nd and Sm iso- 
topes [j-ISl. In this letter, we report a first systematic cal- 
culation of electric dipole {E\) responses for these heavy 
isotopes with the shape phase transition, using a non- 
empirical approach with the Skyrme EDF, namely, the 
quasiparticle-random-phase-approximation based on the 
Hartree-Fock-Bogohnbov ground states (HFB+QRPA). 

GDRs in heavy deformed systems have been previously 
investigated using the separable QRPA with Skyrme 
EDF The separable approximation in the QRPA, 

perhaps, provides a good approximation for the GDR. It 
is, however, difficult to analyze the low-lying states be- 
cause the structure of normal modes is non-trivial and 
significantly affected by the detailed shell structure. The 
low-energy El strengths, which are often discussed as the 
pygmy dipole resonance (PDR) , have attracted a consid- 
erable interests. A nonstatistical distribution of the El 



strengths close to the threshold has a strong impact to the 
astrophysical r-process nucleosynthesis . In addition, 
the PDR is a typical example of exotic collective modes 
expected in neutron-rich nuclei and has been extensively 
studied with the EDF approaches However, the role 
of deformation on the PDR has been studied only for light 
nuclei [l3, 14 1, except for a recent study on Sn isotopes 



with the relativistic EDF [15 1. 



We have developed a new parallelized computer code 
of the HFB+QRPA, which is an extended version of 
that developed in Ref. [16|, to add the residual spin- 
orbit interaction. We expect that the residual Coulomb 
plays only a minor role [13, [l3|- To describe the nu- 
clear deformation and the pairing correlations, simul- 
taneously, in good account of the continuum, we solve 
the HFB equations [18| in the coordinate space using 
cylindrical coordinates r = (p, z, 0) with a mesh size of 
Ap = Az = 0.6 fm and a box boundary condition at 
(Pmax, -Zmax) = (14.7, 14.4) fm. Wc assumc axial and re- 
flection symmetries in the ground state. The quasiparti- 
cle (qp) states are truncated according to the qp energy 
cutoff at Ea < 60 MeV. We introduce the additional 
truncation for the QRPA calculation, in terms of the two- 
quasiparticle (2qp) energy as Ea + Ep < 60 MeV. This 
reduces the number of 2qp states to, for instance, about 
38 000 for the K"" = 0" excitation in ^^''Sm. To calculate 
the QRPA matrix elements and to diagonalize the ma- 
trix, it takes about 390 CPU hours and 135 CPU hours, 
respectively. A detailed description of the method and its 
computational aspects will be soon published elsewhere. 
The similar calculations of the HFB+QRPA for ax ially 
deformed nuclei have been recently reported [H [l9l-El|. 

For a normal (particle-hole) part of the EDF, we em- 
ploy the SkM* functional [22 L For the pairing energy, 
we adopted the one in Ref. 23'] that depends on both 
the isoscalar and isovector densities, in addition to the 
pairing density, with the parameters given in Table HI of 
Ref. [il. 

Since the full self-consistency between the static mean- 
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FIG. 1: Intrinsic electric quadrupole moments of Nd and Sm 
isotopes. For spherical nuclei, values extracted from calcu- 
lated B{E2) are shown by triangles. 



field calculation and the dynamical calculation is slightly 
broken by the neglected residual terms and the trunca- 
tion of the 2qp states, the spurious states may appear at 
finite excitation energies. In the present calculation, the 
excitation energies of the spurious states with K"^ = 0~ 
and 1~, corresponding to the center-of-mass motion, be- 
come imaginary in ^^''Sm, 1.46z MeV and 1.60i MeV, 
respectively. Small contamination of the spurious com- 
ponent does not affect the GDRs because they are far 
apart in energy. The El strength of low-lying dipole 
states might be slightly infiuenced. However, it will not 
affect the conclusion of the present study. 

Figure [T] shows the ground-state deformation of Nd and 
Sm isotopes obtained with the HFB calculation. The cal- 
culated intrinsic electric quadrupole moments are com- 
pared with the experimental values ^23] . The calculation 
well reproduces the evolution of quadrupole deformation 
for N > 86. For spherical nuclei with iV = 82 and 84, we 
also plot the values deduced from the B{E2;0^ 2+) 
obtained by the QRPA calculation. Collectivity of the 
2+ state is apparently overestimated at = 84, because 
these nuclei are so soft with respect to the quadrupole de- 
formation that the QRPA cannot describe the 2+ state 
properly. 

Based on these HFB ground states, we perform the 
QRPA calculation to obtain the excitation energies, hiOi, 
and the transition matrix elements, {i\Flj^\0). The pho- 
toabsorption cross section is calculated as 
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where Ef = {fiiu.f + 7^/4 Q. The smearing width 7 is 
set to 2 MeV, which is supposed to simulate the spread- 
ing effect, F-^, missing in the QRPA. The electric dipole 
operator Pl^^ is defined as Eq. (6-175) of Ref. 0|. 

We show in Fig. [5] the calculated photoabsorption 
cross sections in the GDR energy region together with 
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FIG. 2; (Color online) Photoabsorption cross sections in Nd 
and Sm isotopes as functions of photon energy. The experi- 
mental data 0, Si are denoted by filled squares. 



the available experimental data 0, The GDR 

peak energies well agree with experimental values, and 
produces the deformation splitting in I50'i52]\^jj ^^^^ 
I52,i54gj^ Previously, the separable QRPA calculations 
with SkM* % To'] produce a significantly larger second 
peak for these deformed nuclei, which disagrees with the 
experiment. This is diminished in the present full QRPA 
calculation. The GDR width calculated with 7 = 2 MeV 
is also in good accordance with the experimental values. 
The QRPA accounts for the Landau damping, which is 
a fragmentation of the GDR strength into nearby 2qp 
states, but not for the spreading effect which corresponds 
to a fragmentation into more complex states. The nice 
agreement on the broadening indicates that the smearing 
width 7 = 2 MeV has a good correspondence with the 
spreading width in these nuclei. 

The isotopic dependence of the peak broadening is well 
reproduced, surprisingly, even for the transitional nuclei. 
The width for = 82 and 84 is calculated as F w 4.5 
MeV, and it gradually increases to about 6 MeV for 
N = 88, then the splitting becomes visible for A^ > 90. 
Here, the width F is evaluated by fitting the calculated 
cross section with a Lorentz line. This broadening effect 
is commonly interpreted as the mode-mode coupling ef- 
fects to the low- lying collective modes Q . In the present 
QRPA calculation, the mode coupling is not explicitly 
taken into account. However, the QRPA based on the 
deformed HFB state may implicitly include a part of the 
coupling effect. Figure [H shows that the isotopic depen- 
dence can be well accounted for by the gradual increase 
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FIG. 3: (Color online) Photoabsorption cross sections of 
I'^'Sm and ^^^Sm employing the SkM*, SLy4 and SkP func- 
tionals. The K'^ — 0~ and — 1~ components are shown 
by dotted and thin solid lines, respectively in ^^''Sm. 



of the ground-state deformation. However, a small in- 
crease of the width from ^^^Nd to ^^"'Nd observed in the 
experiment cannot be fully reproduced in the calculation. 
Since the HFB calculation produces the spherical ground 
state for ^*'*Nd, this requires an explicit higher-order cal- 
culation beyond the QRPA. 

We may notice another small disagreement in the peak 
shape: The calculated GDR peak has a shoulder in the 
spherical nuclei, and this shoulder is becoming the third 
peak in the deformed nuclei. This is due to the Lan- 
dau fragmentation, however, this feature is not clearly 
observed in the experiments. As is discussed in the fol- 
lowings, detailed properties of the Landau fragmentation 
depend on the choice of the Skyrme EDF. For instance, 
the fragmentation effect becomes weaker with the SkP 
functional, to give a better agreement with the experi- 
ments. 

Figure |3] shows photoabsorption cross sections in 
I44,i54gj^ obtained by employing the SLy4 [2^ and 
SkP [l^l ■ The experimental GDR peak structure is repro- 
duced not only by using the SkM* functional but also by 
using other commonly-used Skyrme functionals, though 
the SLy4 gives slightly smaller GDR peak energies. The 
SkP reproduces the energy and the shape best among the 
three functionals. For spherical nuclei, the cross section 
obtained with the SkP functional can be nicely fitted 
by a single Lorentzian curve. For deformed nuclei, the 
lower {K = 0) peak shows a Lorentz shape for any of the 
functional, while the upper (K = 1) peak shows visible 
distortion for SkM* and SLy4. This difference may be 
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FIG. 4: Calculated low-energy El strengths B{E1;0+ -5> 1") 
as functions of energy in i'*2,i50,i52]sg-j rpj^^ arrow indicates 
the neutron emission threshold energy. The panel (d) shows 
the unperturbed strengths multiplied by 1/15 in ^^^Nd. The 
neutron threshold of "^Nd is 11.6 MeV. 
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FIG. 5: Same as Fig. [4] but 
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shows the unperturbed strengths multiplied by 1/15 in 
The neutron threshold of "*Sm is 11.9 MeV. 



*Sm. 



due to different properties of the Landau fragmentation 
because the 2qp states in the background are more widely 
spread in energy with a smaller effective mass. Actually, 
the effective mass is largest {m* /m 1) in SkP. 

Next, we discuss the low-energy El strengths. Photon 
scattering experiments at Technische Universitat Darm- 
stadt [2^ reported the 1~ states up to 9.9 MeV and found 
the concentration of the dipole strength in = 82 semi- 
magic nuclei at energies between 5.5 and 8 MeV. Fig- 
ures [3] and [5] show the the El strengths below 10 MeV 
calculated with the SkM* functional. The SLy4 and SkP 
functionals provide very similar results. In the spherical 
^^^Nd and ^''''Sm nuclei, we can see a concentration of 
the dipole strength in between 8 and 10 MeV. The same 
concentration can be seen in the unperturbed strength 
distribution, which may suggest a weak collectivity. Ap- 
parently, the calculated strength distribution is too high 
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FIG. 6: (Color online) Calculated transition densities of the 
K"" = 0" state at 7.49 MeV and of the K'' = 1" state at 7.33 
MeV in ^^*Sm. 



in energy. Similar disagreements have been also observed 
in the relativistic QRPA calculation 



27| . However, the 



B{El t) value summed up to 10 MeV in "''Sm is 0.27 
e^fm^ which reasonably agrees with experimental value 
0.20 e^fm^. It seems to suggest that these low-energy El 
strengths are redistributed to two-phonon (2+ (g)3~) and 
multi-phonon states, which are located at energies 3 ~ 7 
MeV, by higher-order coupling effects [2^. It is quite 
challenging to describe such multi-phonon states on the 
basis of the nuclear EDF method [i^ [2^, however, it is 
beyond the scope of the present analysis. 

In contrast to the spherical nuclei, the deformed nuclei 
show significant fragmentation of the El strength into 
low-energy states. For i50,i52j^|^ ^^^^ ^^^'^^''Sm, the calcu- 
lated energies of the lowest if = 0~ states are 0.97, 1.40, 
1.10, and 1.37 MeV, and those of if^ — 1~ states are 
1.80, 1.93, 1.60, and 1.49 MeV, respectively. These values 
agree well with available experimental values of octupole 
vibrational states js^. The deformation significantly in- 
creases the total low-energy El strength. Actually, the 
summed strengths up to 10 MeV, J2B{E1 t) « 1-5 
e^fm^, are about five times larger than those of spherical 
nuclei. It may be of significant interest to study how this 
low-energy El enhancement due to deformation affects 
the element synthesis scenario. 

We also confirm that the dipole states around 7 MeV 
have a character different from the GDR. In Fig. \6\ the 
transition densities, approximately projected to the lab- 
oratory frame [sij are shown in arbitrary scales. These 
dipole states predominantly have an isoscalar character: 
The proton and neutron transition densities have the 
same sign. The neutron dominance can be seen at r > 7 
fm, however, the proton contribution does not completely 
vanishes. These characters of the FDR are consistent 
with the previous QRPA analysis for spherical nuclei of 
N = 82 isotones 27]. Increasing N/Z ratio away from 
the stability, we expect that the neutron-dominant char- 
acter develops in the surface region. 

In summary, we have investigated effects of the shape 
transition on El strength distribution in the rare- 



earth nuclei, using the newly developed parallelized 
HFB-I-QRPA calculation code with the Skyrme EDF. 
This enables us to simultaneously study both high-energy 
GRs and low-energy collective/non-collective states. The 
typical characteristics of the GDR for shape phase tran- 
sition from spherical to deformed nuclei, especially the 
isotopic dependence of broadening and splitting of the 
GDRs, are extremely well reproduced in the calculation. 
We have also found that the deformation plays a signif- 
icant role for the low-energy El strength distribution: 
The El strength is distributed to low-energy states and 
the total strength at E' < 10 MeV is roughly five times 
enhanced, compared to the spherical nuclei. The low- 
energy strengths in spherical ^^''Sm are calculated too 
high in energy. Inclusion of the higher-order mode-mode 
coupling is desired for further improvements. System- 
atic calculations with the HFB-I-QRPA for spherical-to- 
deformed and light-to-heavy nuclei help us not only to 
understand and to predict new types of collective modes 
of excitation, but also to shed light on the nuclear EDF 
of new generations. 

One of the authors (K.Y.) is supported by the Special 
Postdoctoral Researcher Program of RIKEN. The work is 
supported by Grant-in- Aid for Scientific Research (Nos. 
21340073 and 20105003) and by the Joint Research Pro- 
gram at Center for Computational Sciences, University 
of Tsukuba. The numerical calculations were performed 
on RIKEN Integrated Cluster of Clusters (RICC). 



[1 
[2 

[3; 

[4 
[5] 
[6] 
[7 

[lo; 

[11 
[12 

[13 
[14 
[15 
[16 

[17 
[18 

[19 
[20 
[21 
[22; 



R. G. Parr, Ann. Rev. Phys. Chem. 34, 631 (1983). 

M. Bender et al., Rev. Mod. Phys. 75 (2003) 121. 

M. V. Stoitsov et al, Phys. Rev. C 68, 054312 (2003). 

A. Bohr and B. R. Motteleson, Nuclear Structure, vol. II 
(Benjamin, 1975; World Scientific, 1998). 

M. N. Harakeh and A. van der Wounde, Giant Reso- 
nances (Oxford, 2001). 

B. L. Berman and S. C. Fuhz, Rev. Mod. Phys. 47, 713 
(1975). 

P. Carlos et al, Nucl. Phys. A172, 437 (1971). 

P. Carlos et al, Nucl. Phys. A225, 171 (1974). 

V. O. Nesterenko et al, Phys. Rev. C 74, 064306 (2006). 

V. O. Nesterenko et al, Int. J. Mod. Phys. E 17, 89 

(2008). 

S. Goriely and E. Khan, Nucl. Phys. A706, 217 (2002). 
Paar et al, Rep. Prog. Phys. 70, 691 (2007). 
Yoshida, Phys. Rev. C 80, 044324 (2009). 
Ebata et al, arXiv: 1007.0785 

P. Arteaga et al., Phys. Rev. C 79, 034311 (2009). 
Yoshida and N. V. Giai, Phys. Rev. C 78, 064316 



(2008). 

J. Terasaki et al, Phys. Rev. C 71, 034310 (2005). 
J. Dobaczewski et al, Nucl. Phys. A422, 103 (1984). 
S. Peru and H. Goutte, Phys. Rev. C 77, 044313 (2008). 
C. Losa et al., Phys. Rev. C 81, 064307 (2010). 
J. Terasaki and J. Engel, arXiv: 1006.0010 
J. Bartel et al., Nucl. Phys. A386, 79 (1982). 



5 



[23] M. Yamagami et al, Phys. Rev. C 80, 064301 (2009). 
[24] http://ie.lb l.gov/ensdf/ 

[25] E. Chabanat et al, Nucl. Phys. A627, 710 (1997). 
[26] A. Zilges et al, Phys. Lett. B 542, 43 (2002). 
[27] N. Paar et al, Phys. Rev. C 67, 034312 (2003). 
[28] D. Sarchi et al., Phys. Lett. B 601, 27 (2004). 



[29] E. Litvinova et al, Phys. Rev. Lett. 105, 022502 (2010). 
[30] P. C. Sood, D. M. Headly and R. K. Sheline, At. Data 

Nucl. Data Tables 47, 89 (1991). 
[31] D. Pena Arteaga and P. Ring, Phys. Rev. C 77, 034317 

(2008). 



